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ABSTRACT 

We present new integral field spectroscopy of the gravitationally lensed broad ab¬ 
sorption line (BAL) quasar H1413-I-117, covering the ultraviolet to visible rest-frame 
spectral range. We observe strong microlensing signatures in lensed image D, and we 
use this microlensing to simultaneously constrain both the broad emission and broad 
absorption line gas. By modeling the lens system over the range of probable lensing 
galaxy redshifts and using on a new argument based on the wavelength-independence 
of the broad line lensing magnifications, we determine that there is no significant broad 
line emission from smaller than ~20 light days. We also perform spectral decomposi¬ 
tion to derive the intrinsic broad emission line (BEL) and continuum spectrum, subject 
to BAL absorption. We also reconstruct the intrinsic BAL absorption profile, whose 
features allow us to constrain outflow kinematics in the context of a disk-wind model. 

We find a very sharp, blueshifted onset of absorption of 1,500 km/s in both Civ and 
Nv that may correspond to an inner edge of a disk-wind’s radial outflow. The lower 
ionization Siiv and Aim have higher-velocity absorption onsets, consistent with a de¬ 
creasing ionization parameter with radius in an accelerating outflow. There is evidence 
of strong absorption in the BEL component which indicates a high covering factor for 
absorption over two orders of magnitude in outflow radius. 
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Introduction 


The central engines of powerful Type I active galactic nuclei (AGNs) contain gas that is moving 
at extreme velocities, apparent from their broad emission lines (BELs) and, in the case of broad ab¬ 
sorption line (BAL) quasars, in blueshifted absorption troughs. These spectral features arise from 
compact gas flows close to the accretion dis k, apparent from the pres e nce of high ionization lines 


and s ize scales from reverberation mapping (iBlandford fc McKeelll982l : IPetersonlll99,‘ll : iKaspi et al 


20071 1. These compact flows are probably the origin of larger-scale outflows (e.g. Chelouche & Net- 


zer 2005; Crenshaw &: Kraemer 2007; Storchi-Bergmann 2010) and result in feedback that is believed 
to be important on many scales, regulating supermassive black hole (SMBH) fueling, quenching star 
formation, and so driving galaxy-SMBH coevolution (e.g. Silk and Rees 1998; Fabian 1999; Som- 
merville et al. 2008; Zubovas &: King 2012) as well as influencing the intergalactic and intracluster 
medium and affecting galaxy and large-scale structure formation (e.g. Scannapieco & Oh 2004; 
Rafferty, McNamara &: Nulsen 2006; Sommerville et al. 2008; Fabjan et al. 2010; Wang 2015). Yet 
despite the importance of these processes, the connection between compact and large-scale outflows 
is not understood. In the case of the broad emission line region (BELR), even the basic kinematics 


(Kasoi et al. 

2000. 

2097; 

Dennev et al. 

2009: 

Bentz et al. 

2009: 


Joroshenko et al 


crole nsing results have been consistent with both viria l BELBs (lO’Dowd et al 


20 111 ) and with more complex, multi-component flows (jBichards et al 


2004 


20 


2011 


2), and mi- 


Guerras et al 


Shise et al. 

201 oi l. 


BAL quasars offer perhaps the clearest evidence of nuclear outflows in AGNs. Approxi- 
mately 15% of optica lly selected quasars exhibit high velocity, blueshifted absorption features 


(jReichard et al.l 1200,01 ). The outflows responsible for these features may be common to most 


quasars if this 15% represents a special liii e of sight to outflows with a similar covering factor 


(Wevmann et al 


1991 


Gallagher et al.ll2007l ). or a particular phase in the lifecycle of quasars. An 


accurate measurement of the absorption profile of a BAL quasar would go a long way to revealing 
the kinematics of these outflows, however this measurement is complicated by the fact that their 
intrinsic absorption profile cannot be recovered without knowing the shape of the underlying BEL 
that it absorbs. These BEL profiles can vary dramatically between quasars, and in the case of BAL 
quasars may be strongly affected by a biased viewing angle. 

The gravitationally lensed BAL quasar H1413-I-117 (the Cloverleaf; Magain et al. 1988) offers 
a unique opportunity to simultaneously study both BEL and BAL flows in the same source. The 
four lensed images of this z = 2.55 quasar exhibit deep BAL features in high-ionization lines, as 
well as low-ionization absorption in Aim, making it a LoBAL quasar. Signif i cant microlensing has 


Monier et al. 

1998; 

Aneuita et al. 

2008: 

Hutsemekers et al. 


201 




In microlensing, substructure 
within the foreground lensing galaxy causes changes in the magnification of individual lensed images. 
The degree of microlensing depends on size of lensed region, and in the case of H1413-I-117 this 
has resulted enhanced magnification in the continuum in lensed image D compared to the broad 
emission lines in this component. 
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Given a model of the microlensing magnification of the lens system, these differential magnifi¬ 
cations can be used to derive size constraints of emission regions. This approach has been used to 
contrain accretion disk size and temperature gradient in several lensed quasars (e.g. Kochanek et 
al. 2004; Bate et al. 2008; Floyd, Bate Sz Webster 2009; Blackbu rne et al. 2011; Jimenez- Vicente 


et al. 2014; Blackburne et al. 2015), including in H1413-I-11 7 by iHutsemekers et al.l feOlOjl 


after HIO). It has also provided constraints on BELR size (IKochaneB l2004l : 


Wavt 

1 . O’Dowd & Webstei]l2005 

Sluse et a . 

to 

0 

0 

O’Dowd et al. 

2011: Sluse et al. 

2011: 


2007; 


Guerras et al 


.II2OI 3I 1 


here- 


Richards et al 


2004: 


and kinematics (Richards et al 


Guerras et al.ll2013l l in a number of quasars, however 


no such constraints have been derived for the BELR of H1413-I-117. 


Another novel use of microlensing takes advantage of differential magnifications between lensed 
images to algebraically deco m pose the spect r al con tributions of different physical regions within 
the quasar (jSluse et al.ll2007l ). lAngonin et al.l (jl990|) (hereafter ADO) and HIO apply this technique 
to H1413-I-117, shedding light on the true shape of the broad absorption profile, and in the case 
of HIO decomposing into microlensed and non-microlensed components which reveal substructure 
in each component and may point to a two-component outflow comprising a polar outflow and 
equatorial disk-wind. 


In this paper we present new Gemini North GMOS integral field spectroscopy of H1413-I-117, 
and based on the strong microlensing observed in lensed image D we derive new statistical size 
constraints on the BELR. We also present a new spectral decomposition, deriving BELR and 
continuum components, as well as the intrinsic BAL profiles for Nv, Siiv, Giv, and Aim. 

In Section [2] we describe the observations and extraction of lensed image spectra. In Section 
[3] we analyse the microlensing signatures in this new epoch of data. In Section [3.31 we use lensing 
simulations to explore the size dependence of the observed microlensing signatures and deriving size 
constraints for the BELR. In Section 0] we decompose the spectra in terms of the broad emission 
line region and continuum, emphasizing the physical interpretation of the decomposed regions, and 
we derive the intrinsic BAL absorption profile. We discuss the interpretation of these decomposed 
components in the context of disk-wind outflow models in Section [5l In Section [6] we summarize 
our conclusions. 


2. Observations and Spectral Extraction 


H 1413-H17 was observed on the night of April 5th, 2006 with the GMOS Integral Eield Unit 
(lEU) (|Allingtou-Smith et al.ll2002l : iHook et al.ll2004l ) on the Gemini North telescope (program ID: 
GN-2006A-Q-15). 3x20 minute exposures were taken in one-slit mode using the B600 grating with 
a spectral range of 396lA to 6845A and a resolution of i? = 1688. At the redshift of H1413-I-117, 
2 : = 2.55, this covers the Lya through Giii] broad emission lines. Seeing conditions were good, with 


a PSE EWHM of 0'.'8. 
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Fig. 1.— 2-D image created from collapsed IFU data cube. Letters indicate the standard designa¬ 
tion of quasar lensed images in H1413-I-117. 


2.1. Data Reduction 


The data was reduced using the Gemini IRAF package vl.lO under IRAF 2.14.1. We used 
calibration files taken as part of the standard Gemini Facility Calibration Unit (GCAL) suite. 
Mostly, the standard CMOS IFU reduction procedure was followed for overscan subtraction and 
trimming, bias subtraction, flat fielding and response correction, arc calibration, fiber extraction, 
sky subtraction, flux calibration, cube reconstruction with atmospheric dispersion correction. 


Additionally, we corrected the GCAL flat and the twilight flat for scattered light by fitting 
a surface to residual light detected between the fiber bundles. That correction was not necessary 
for the science observations, the signal being too weak to generate any significant scat tered light. 
The cosmic rays were corrected before extraction with LA-Cosmic (|van Dokkumll200lh . followed, 
after extraction, by a floor and ceiling masking of remaining extreme hits. The calibrated cubes for 
each science observation were aligned and combined with a customized version of the pyfu software 
package^. 


^Scattered light correction tasks by Bryan Miller, Gemini Observatory (version 2010). LA-Cosmic IRAF wrapper 
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Figured] shows the reduced data cube collapsed into a 2-D image. All lensed components are 
distinct and well-resolved. 


2.2. Extraction of Image Spectra 


Lya N V 01 C11 Si IV CIV He II AI III C III] 



Fig. 2.— Extracted spectra for the four quasar lensed components: black: image A, green: image 
B, blue: image C, red: image D. Vertical dotted lines mark the positions of emission lines assuming 
z = 2.55. Dashed lines indicate the best linear fits to the continuum spectra of lensed image A & 

D. 


The lensed images are well resolved in these observations, allowin g for a good extracti on of 
their spectra. We perform this extraction using the method described in lo ’Dowd et al. ( 2011 ). We 
refer the reader to that paper for a detailed description, however in short: relative lensed image 
fluxes at each wavelength bin are determined by least-squares fitting of reconstructed Point Spread 
Functions (PSFs). A separate PSF is determined for each lOOA-wide spectral bin. Each given PSF 
is initially estimated using carefully-masked fragments of the IFU frame (stacked over lOOA), and 


and pyfu suite by James Turner, Gemini Observatory (version 2010). All, except the pyfu software, are now part of 
the public Gemini IRAF package (vl.13). The pyfu software is available on the Gemini Data Reduction Forum at 
http: //drforum.gemini.edu/topic/pyfu-datacube-mosaicking-package / 
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are iteratively improved with a series of least-squares minimizations of each frame. 

Figure [2] shows the extracted spectra of the four quasar images. We see distinct broad emission 
lines for Lya, Nv, Oi, Cii, Siiv, Civ, and Cm] as well as prominent iron emisison features at ~1600- 
I 8 OOA. Emission from Hell and Aim is also apparent, blended with iron and Cm] respectively. 
Clear redshifted broad absorption troughs are seen for Lya/Nv, Siiv, Civ, and Aim. 


3. Microlensing Analysis 

Gravitational leasing is achromatic, and so in the absence of microlensing, differential extinc¬ 
tion, or time-dependent spectral variation, the spectra of leased images of a multiply-imaged quasar 
are identical modulo a scaling factor. Microlensing is typified by a difference in the magnification 
of one or more leased images compared to the time-averaged magnification for that leased image 
(which we will refer to as the macrolens magnification, or just macro-magnification). Differential 
microlensing is characterized by magnifications that vary between different emission regions, result¬ 
ing in a change in the shape of the spectrum of the affected leased image. Differential microlensing 
is sensitive to emission region size and projected location, is frequently seen between BELs and 
continua, and also often shows wavelength dependence within these features. 


3.1. Intrinsic Variability and Dnst Extinction 


Spectral variability coupled with different light travel times can result in spectral differences 
between leased images. In the case of H1413-I-117 the time del ay is potentially significant, es - 
timated, for example at 23.4±4 days between images A and D (iGoicoechea fc ShalvapinI 120101 ). 
This variability can potentially lead to differences in four key observables between leased images: 
BAL and BEL profile shapes, continuum-to-BEL flux ratio, and the spectral index of the contin¬ 
uum. While it is difficult to extricate the differences due to intrinsic variability from those due to 
microlensing, we can comment on the effect of this on our analysis. 


Variation in broad absorption and emission line profile shapes would be very difficult to dis¬ 
criminate from microlensing. Eortunately the effect is probably minimal. As we show in Section [32] 
the BEL profiles are identical between leased images, precluding both differential microlensing and 
intrinsic variability in the BELR spectrum. It is not so easy to compare BAL profiles, however we 
note that, while this profile has been observed to vary, this variation occurs over months and years, 
not days (HIO). As we will see in Section 13.31 we estimate the size of the BELR (and hence also the 
BAL wind) at >100 light days, and so their intrinsic variation timescales should be significantly 
longer than the time delay. 


Variation in the continuum-to-BEL ratio between lensed images cannot be ruled out, and we 
do see significant variation in this ratio, which we attribute to microlensing because it is consistent 
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with the microlensing behavior of the lens in recent years. In Section [J] we use this difference in 
continuum-to-BEL ratios to decompose the spectrum into intrinsic components. It’s important to 
note that for the purpose of the decomposition of spectral components, it does not matter why this 
ratio varies between images—whether by microlensing or intrinsic variability—the decomposition 
is still valid as long as the BELs don’t show wavelength-dependent variations. 

Variation in the slope of the continuum could also be due to intrinsic variability, and we do 
indeed observe a difference in slope (Sect. 13.21) . however we don’t investigate this in detail in this 
paper. 

Differential extinction between lensed images will not affect the continuum-to-BEL ratios as 
our measurement of the continuum strength is made at the relevant BEL wavelength range. Nor 
will it significantly affect broad line profiles, which appear identical between lensed images anyway. 
Differential extinction only affects the continuum slope and ratios between BELs, which are not a 
focus of his work. 


3.2. Flux Ratios and Microlensing Wavelength Dependence 

To inspect differential microlensing properties of H1413-I-117 we first obtain an estimate of the 
unabsorbed continuum spectrum for each lensed image by performing a linear fit to regions of the 
spectrum least affected by emission or absorption (Fig. [2] shows fits for images A & D). Prominent 
iron emission features force us to choose small windows at rest frame 1420—1460A, 1685—1692A, 
and 1795—1805A. While the continuum is unlikely to be linear over the full range of this spectrum, 
we find that the fit gives a very plausible estimate, particularly redward of the Nv line. Subtraction 
of the fitted continua allow us to measure the average flux ratios across each BEL for each lensed 
image with respect to image A. These ratios and the average flux ratios for the continuum at each 
line center are given in Table [TJ For the purpose of inspection of the microlensing properties, we 
find that changing the continuum level and slope by allowing variation of ±5% (the width of the 
pale gray region surrounding fits in Fig. [2|) does not significantly affect results. 

Figure [3] (upper panels) compares the spectra of all lensed components around the broad emis¬ 
sion lines. To highlight the effect of differential microlensing we scale the spectra of lensed images 
B, C, & D so that the average continuum emission at each line (from the continuum fits) matches 
that of image A. We see significant differential microlensing between broad emission lines and con¬ 
tinuum in lensed image D versus all other images, with the emission lines at lower magnification 
in D relative to its continuum. Other lensed images show weak differential magnification between 
lines and continuum. Due to an unknown lens redshift the lens model is not well established and so 
the macro-magnifications are not well known. For this reason we don’t know which, if any of A, B, 
or C are free of microlensing. Nonetheless, we can conclude that image D is undergoing significant 
differential microlensing between BELR and accretion disk. 

These results are consistent with observations over the past decade (Angonin et al. 1990, 
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Fig. 3.— Upper: Spectra at broad emission lines with fluxes scaled so that the continua match in 
regions bounding lines - black: image A, green: image B, blue: image C, red: image D. Lower. 
Ratios of the spectra of lensed images A &: D at broad emission lines after subtraction of best-fit 
continuum spectra (black). Blue shows the continuum-subtracted broad lines for lensed image A. 
Green lines in the lower panels show the average broad-line flux ratio for each line (Tabled]). Orange 
lines indicate the BEL rest-frame wavelengths. 


Ostensen et ah 1997, Monier et al. 1998, Anguita et al. 2008, HIO), which show strong anomalous 
line-continuum flux ratios in image D compared to A, B, and C. 

Our observations show wavelength-dependent differential microlensing of the continuum in the 
D/A and C/A ratios. The D/A ratio tends to more extreme values (deviating from the average 
BEL ratio) towards shorter continuum wavelength, as was also observed by HIO. It is tempting to 
interpret this as a simple inverse correlation between emission region size and microlensing intensity, 
however we note that the observed C/A ratios have the opposite behaviour, tending towards the 
average C/A BEL ratio towards /on^er wavelengths. In Section 13.31 we show that we cannot expect 
a monatonic trend of magnification ratio with size across multiple orders of magnitude in emission 
region size. Proper analysis of differential microlensing across the continuum requires careful lensing 
simulation, however we do not perform that analysis in this work. 

To investigate the microlensing state of the BELs we determine the wavelength-dependent flux 
ratios for the BEL SEDs for lensed images A and D. These are just the ratios of the spectra of 
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B/A C/A D/A 


line continuum line continuum line continuum 


Nv 

0.78 

± 

0.02 

0.92 

± 

0.02 

0.93 

± 

0.03 

0.81 

± 

0.02 

0.45 

± 

0.03 

0.77 

± 

0.03 

Silv line 

0.82 

± 

0.02 

0.91 

± 

0.01 

0.93 

± 

0.02 

0.79 

± 

0.01 

0.38 

± 

0.02 

0.75 

± 

0.02 

Civ 

0.81 

± 

0.01 

0.91 

± 

0.01 

0.85 

± 

0.01 

0.77 

± 

0.01 

0.41 

± 

0.01 

0.73 

± 

0.01 

Cm] 

0.84 

± 

0.02 

0.88 

± 

0.02 

0.87 

± 

0.03 

0.71 

± 

0.02 

0.38 

± 

0.03 

0.68 

± 

0.02 

Fe complex 



0.79 ± 0.02 





0.84 ± 0.01 





0.41 ± 0.01 




mean BEL 

0.80 ± 0.04 

0.89 ± 0.04 

0.40 ± 0.02 

mean cont. 

0.91 ± 0.03 

0.77 ± 0.02 

0.73 ± 0.02 


Table 1: Measured flux ratios for broad emission lines and continuum at BEL locations for lensed 
images B, C, and D with respect to A, based on linear fit to continuum (see Sect. [ 321 ), and avoiding 
absorption features. Also included are mean BEL and continuum flux ratios. 


images A and D after subtracting the continuum estimates described above. In Figure [3] (lower 
panels) we can see that the flux ratios for the Nv, SiiV, CiV, and Cm] broad emission lines are 
flat over the wavelength range redward of BAL absorption. Lya is omitted due to its overlapping 
the Nv absorption trough, which makes continuum subtraction impossible. This flatness indicates 
that there is no differential microlensing across the BELR velocity structure redward of the BAL 
absorption, and very little difference between these emission species. The mean and individual BEL 
flux ratios are given in Table HI Error intervals include standard error of the mean and uncertainty in 
the continuum fits. Variation of the level of continuum subtraction up to its estimated uncertainty 
at the location of each emission line does not change the observed flatness of these flux ratios. 


3.3. Comparison to Q2237+0305 and Size Constaints on the BELR 


In lO’Dowd et al.l (l201ll l we analyzed microlensing signatures in the Civ and Cm] BELs of 
Q2237+0305, the Einstein Cross. In that case there was clear differential microlensing across the 
velocity structure of these broad lines, with magnification ratios varying by ~25-30% between line 
center and wings. The BEL magnification ratio approached the highly anomalous continuum mag¬ 
nification ratio at the highest velocities. Our simple model indicated that this trend was consistent 
with the broad line flow being a Keplerian field, and not supporting a pure radiation-driven out¬ 
flow without a substantial gravitational component. This differential microlensing signature sharply 
contrasts the wavelength-independent BEL magnification ratios seen in H1413-I-117 (Sect. [3]). This 
wavelength independence suggests that the BELR in H1413-I-117 is large enough to average over 
fluctuations in the magnification surface. Given the very low redshift of the lensing galaxy in 
Q2237-I-0305 {zl = 0.04) and its commensurately large fluctuation scale-length, differential mi¬ 
crolensing of its BELR is not surprising. To explore the implications of this difference in more 
detail we compare the expected microlensing behavior of these two quasars as a function of source 


size. 
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Although the lensing galaxy of H1413+117 does not have a conhrmed redshift, it is possible 
to explore plausible microlensing scenarios usi ng estimates from literature: zj, = 0.94 based on the 
probability distribution from lens parameters (jMosauera Kochanekll201lh and supported by th e 
zl ~ 0.9 pho tometric redshift of overdensities in the hcM (IKneib et al.lll998l : iFaure et al.l 120041 '): 


and zl = 1.4 ( Surdei et al. 1993 : Witt. Mao k, Schechtei 1995| b tha t is also comparable t o the z 


1.5 absorption systems in the quasar spectrum (jMagain et al.lll988l : lAngonin et al.lll99flli. For the 
H1413 +117 lens macro-model we use the convergence and she m estimates oflWitt. Mao fc Schechter 
(1995) (/? = 1 case), which are nearly identical to those of Mediavilla et al. ( 2009l b We assume 
an 80% smooth matter percentage, consistent with the values estimated in other lenses where 


the line of sight does not pass through the galactic core (Keeton et al. 2006: Chartas et al. 2009: 

Poolev et al. 

o 

o 

Dai et al. 

o 

o 

Bate et al. 

2011 

). Eor Q2237-I-0305 we use the macro-model of 

Trott et al. ( 

201CL and Trott &: Wavth fpriv. comm.L In this case lensed images are clustered 


around the star-dominated g alactic core and microlensing is consistent with a 0% smooth matter 
percentage (jBate et al.ll201in . which we employ here. We obtained simu lated maps of microlensing 
magnifications from the GERLUMPH database (jVernardos et al.ll2014l l based on these models. 


We simulate the expected microlensing behavior of these lenses, using Monte Carlo simulations. 
Emission region surface brightness distributions are simulated with Gaussian profiles whose half- 
light radii represent emission region sizes. Note that except in the case of very strong microlensing 
events, half -light radius rather than the detailed surface brightness distribution governs microlensing 
properties (jMortonson, Schecter &: Wambseanssl l2005l b A pair of random locations on Gaussian- 
convolved magnification maps (one map for each of a pair of lensed images) yields a possible 
magnification ratio between those lensed images. By varying half-light radius and location we can 
then explore microlensing behavior as a function of emission region size. We determine D/A ratios 
for H1413-I-117 at 50,000 location pairs and B/A ratios for Q2237-I-0305 at 25,000 location pairs. 
We explore sizes from 3 x 10^^ m to 2 x 10^® m, constrained by the resolution and size of the 
GERLUMPH maps. 


Eigure [4] (top 3 panels) shows magnification ratio versus emission region size for a random 
selection of 26 of these locations. Magnification ratios are normalized to unity at the macro¬ 
model magnification ratios (the large-scale average ratios). These figures, which we will refer to as 
convergence plots, illustrate two important points: 


1 ) As expected, magnification ratios will converge on the macro-model ratios above a certain 
emission region size (the convergence size). We see that for H1413-I-117 (top 2 panels) the D/A 
convergence size is ~ 5 x m (~20 light days) for both the zl = 0.94 and zl = 1.4 cases. If the 
wavelength independence that we observe in the BEL magnification ratios (Sect. 13.2h is due to such 
convergence, then this gives us an approximate lower limit on the size of observed BELRs. In the 
case of Q2237-I-0305 (panel 3) the convergence of the B/A magnification ratio does not occur until 
significantly greater than 10^^ m, beyond the limit of our simulated magnification maps, and so the 
observation of differential microlensing across the full velocity width of its BELR is not surprising. 
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Fig. 4.— Top panel: Simulated magnification ratios as a function of source size for H1413+117, 
lensed images A and D, with an assumed lens redshift zl = 0.94 (showing 26 of 50,000 simulated 
tracks). 3^^ panel: As above, but for zl = 1-4. 3^^ panel: As above, but for Q2237+0305 

{zl = 0.04; showing 26 of 25,000 simulated tracks). Bottom panel: Percentage of simulated tracks 
that are flat as a function of maximum emission region size. Flatness is defined as < 5% variation 
in magnification ratio over 0.5 dex in size up to the maximum emission region size. The red dotted 
line shows the half-light radius of a Shakura-Sunyaev thin disk accreting at the Eddington limit 
around a non-rotating black hole with mass 4 x 10^ Mq (left line) and 10® Mq (right line), observed 
at 150 nm. The red dashed line is the 150 nm continuum sourse size estimate of HIO. 
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2) We expect tracks on the convergence plot to be flat for sizes larger than the convergence 
size, however we also see many flat tracks for very small emission regions. This independence of 
magnification ratio with size often occurs in cases where the magnification ratio deviates signifi¬ 
cantly from the macro-model value (i.e. where there is strong microlensing). For H1413-I-117, most 
simulated tracks are relatively flat below ~ 10^^ m. There is an intermediate range over which all 
tracks converge on the macro-model ratios—the convergence range. The convergence range spans 
10^^ m to ~ 5 X 10^^ m for H1413-I-117, and in it we see strong size dependence of magnification 
ratios, and here the only flat tracks are those already close to the macro-model ratio. 

The interpretation of this behavior is straightforward: The small-scale flatness seen the conver¬ 
gence plots occurs when a small emission region falls in relatively featureless regions between strong 
magnification gradients in the caustic network. As size increases into the convergence range, our 
larger source is beginning to overlap high or low magnification structures in the caustic network, 
resulting in size-dependent, differential microlensing. When the emission region is large enough to 
sample many such structures, tracks again flatten at the macro-model magnification ratio. 

The magnification ratios of the Civ and Cm] BELs in H1413-I-117 vary by less than 5% across 
the unabsorbed component of their velocity structure (Sect. 13.2p . and so must reside on a “flat” 
region of a convergence plot track. We can place more formal constraints on the possible sizes of 
these emission regions by determining the number of tracks that show similar flatness over a given 
size range. We conservatively assume that the broad line gas flow spans a size range of 0.5 dex 
across its entire velocity structure. In fact the BELR probably spans a larger range, but choosing 
a smaller range yields more tracks within our flatness constraint, and so produces weaker (more 
conservative) size constraints. 

Eigure [4] (bottom) shows the percentage of tracks consistent with <5% variation in the magni¬ 
fication ratio over 0.5 dex in emission region size as a function of the maximum size of the emission 
region, for both H1413-I-117 {zl = 0.94) and Q2237-I-0305. 

In the case of Q2237-I-0305 Figure [4] shows a reasonable likelihood of observing size-independent 
microlensing over this full range of emission region sizes. For H1413-I-117, size-independent mi¬ 
crolensing is unlikely (< 5%) for BELRs spanning 0.5 dex with maximum size in the range ~ 10^^ m 
to ~ 10^^ m (< 1% for ~ 2 X 10^^ m to 5 X 10^^ m). At larger sizes, size independence will become 
more probable, and extrapolating our results we estimate that size independence will become likely 
(>50%) at BELR sizes of 2-3x10^^ m, or around 100 light days. It is important to note that HIO 
also observe this size-independent microlensing of the BELR at a different epoch, greatly decreasing 
the probability that the BELR is in the 10^^ m to 10^® m range. For a reasonable probability 
of observing size independent magnification ratios in the BELR of H1413-I-117 in two independent 
epochs, the BELR probably needs to be larger than 10^^ m (~40 light days), and any significantly 
luminous component likely arises well above the convergence size of ~ 5 x 10^^ m (~20 light days). 
These constraints are for the zl = 0.94 case; larger values of zl push us to larger BELR sizes. 

We can also consider the possibility that the entire BELR is emitting at sizes smaller than 
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I'N-' 10^^ m ( I'N-' 4 light days). We think this unlikely for two related reasons: 1) our observations 
and those of HIO show wavelength-dependent differential microlensing in the continuum D/A ratio, 
which precludes flat tracks at small size needed for size independence in the BEL magnihcation 
ratios; 2) for a high likelihood of size independence at small scales at two epochs, the emission 
region must be closer to a few times 10^^ m, which is approaching the scale of the accretion 
disk (see Sect. 1^. and yet the D/A ratios indicate very strong differential microlensing between 
continuum and BELs at multiple epochs (see Sect. [ 321 ). 


F rom the reverberation mapping BELR size-luminosity relation for luminous quasars of lKasoi et al 
( 2007 1. H1413-I-117 (with its A Fa(1350A)= 4 x 10"^^ erg s“^) is expected to have a Civ BELR in 
the range of 25-90 light day s . From their microlensing study of of 16 lensed quasars (not including 
H1413-I-117), I Guerras et al.l ()2013l ) hnd a BELR size range for high-ionization lines in their lumi¬ 
nous subsample of 36^^4 light days. Both of these results support our estimate of a BELR that is 
signhcantly larger than 20 light days, and very unlikely to be <4 light days. 


3.4. Black Hole Mass and Accretion Disk Size 


We e stimate the black hole mass in H1413-I-117 at 4 x 10^ Mq based on the Civ width and 
using the IVestergaard fc Peterson! (2006j) Ms^f-FWHM relation. We measure the Civ FWHM of 
3,400 km/s based on the unabsorbed, red side. Note that the intrinsic Civ BEL is either symmetric 
or shows a blue excess, as discussed in Section 15.21 and so the correct FWHM and corresponding 
black hole mass may be a larger. 

In order to obtain an upper limi t estimate on the size of th e accretion disk in H1413-I-117 we 
assume a Shakura-Sunyaev thin disk ( Shakura Sz Sunvaevlll973l l with a 10® Mq non-rotating black 
hole accreting at the Eddington limit, and observed at 150 nm. This yields a source half-light radius 
of 4 X 10^^ m, shown in Figure 12 For our black hole mass estimate of 4 x 10® Mq and the same 
accretion parameters the 150 nm source size is 2.7 x 10^® m. Based on differential microlensing 
across the continuum, HIO estimate a continuum source size of 3.4 x 10^® m. 


The accretion disk for both quasars is much smaller than the convergence range. For H1413-I-113 
both size-dependent and size-independent magnihcation ratios are reasonably common at small 
radii, and it’s not surprising that the continuum shows some differential magnihcation in both our 
observations and those of HIO. This is also true for Q2237-I-0305, however size dependence occurs 
less frequently at these small size scales, in ~20% of cases. 

The D/A magnihcation ratio in the H1413-I-117 continuum tends towards the average D/A BEL 
ratio towards longer continuum wavelengths (although never getting close to the BEL ratio), and 
so it’s tempting to interpret this as the accretion disk becoming large enough to start averaging 
over causic structures and hence approach the macro-model magnihcations. However Figure [2 
shows that even a large accretion disk is much smaller than the convergence range. In this regime, 
magnhcation ratio can have a size dependence that either approaches the macro-model value or 
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moves away from the macro-model value, even if all tracks ultimately converge on the same value. 
In short, we cannot interpret a trend towards the macro-model magnifications as a reliable proxy 
for increasing emission region size except in the convergence range. It is therefore not surprising 
that the C/A continuum ratios tend in the opposite direction to the D/A ratios, with shorter 
wavelengths approaching the C/A BEL ratio. 

Proper interpretation of the differential microlensing in the continuum requires higher resolu¬ 
tion magnification maps and a careful statistical analysis, which we leave to a later study. 


4. Algebraic Decomposition of Intrinsic Spectral Components 
4.1. Decomposition of Components 


Lyg NV 0 I C II Si IV 


C IV He II 


Aim CIIIl 
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Fig. 5.— Upper. Absorbed continuum spectrum of lensed image D in H1413-I-117 {Cd from Eqn. [3l 
black). Green dashed line shows the continuum fit to image D. Lower, absorbed broad emission line 
spectrum of lensed image T) {jCd from Eqn. 01 black). The green line shows the observed spectrum 
of image D with cont i nuum fit subtracted. The blue line shows the iron emission template of 
Vestergaard &: WilkesI (|200in with arbitrary scaling. 


We can consider the spectra of the lensed images in H1413-I-117 to be linear combinations 
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of intrinsic spectral components from different emission regions, each subject to absorption and 
lensing that may vary between lensed images and emission regions. If the relative strengths of the 
emission components vary between one lensed image and another due to differential microlensing or 
intrinsic variability, then with some assumptions it can be possible to solve for these components. 


Angonin et al.l (Il990l ) first apply this approach to H1413+117 to derive an intrinsic BAL ab¬ 
sorption profile for Siiv, Civ,and Aim. They describe the spectrum as composed of physical 
co mponents: a cont inuum -|- BEL. HIO explore the approach in more detail, applying the method 
of Sluse et al. ( 2007 1 to H1413-I-117 to extract Lya-I-Nv, Siiv, Civ BAL troughs. In the method, 
the quasar spectrum is explicitly described in terms of a component that is only macrolensed and 
one that is both micro- and macrolensed. Both approaches lead to similar results, and this is be¬ 
cause their assumptions are similar: that the continuum is microlensed but the BELR is not; as a 
result the micro-|-macrolensed component of HIO is essentially the absorbed, lensed continuum. 


Here we derive a new spectral decomposition of H1413-I-117. Our approach is mathematically 
similar to both Angonin et al. and HIO, and we discuss their connection in 14.21 We choose to 
formulate the decomposition in terms of physical emission regions rather than by their degree of 
microlensing. We do this because: 1) our observation of almost constant magnification ratios across 
the emission line velocity structures (see Sect. [3|) ensures a straightforward extraction of the BEL 
component, and 2) we include in our formulation a separate BAL absorption term for both the 
continuum and BELR, which we feel adds clarity to the interpretation of the derived components, 
and 3) any source of variation—microlensing or intrinsic variability coupled with time delay— 
enables a decomposition, and yet the possible presence of the latter means that results can not be 
easily interpreted as microlensed/non-microlensed components. 

We describe the spectrum Ex of lensed image X as a linear combination of two components: 
1) Cx'- the intrinsic continuum spectrum Fc multiplied by the BAL absorption profile affecting 
the continuum Ac and the (potentially wavelength-dependent) magnification mc^x affecting the 
continuum in image X, and 2) Cx- a broad emission line spectrum Fl with its own, potentially 
different absorption profile Al and lensing magnification ruL^x- The F and A components are 
intrinsic to the quasar and identical for all lensed images (modulo quasar variability, but see below). 
The lines of sight to different lensed images do not diverge significantly within the BAL wind (see 
Sect. 15.41) and so all lensed images see the same absorption from this outflow. However separate 
absorption terms are required becaues the accretion disk and BELR are viewed through different 
regions of the BAL wind. The magnification factor is different for each lensed image and in the 
presence of microlensing this factor may differ between the continuum and BELs, and each of these 
may have wavelength dependence. We have neglected to include an explicit term for reddening or 
to account for intrinsic variability in the quasar. Both of these effects cause wavelength-dependent 
multiplicative changes to component spectra, and so in terms of the decomposition we can consider 
these to be folded into the magnification terms; the origin of these multiplicative effects does not 
affect the decomposition. Eor lensed images 1 and 2 we have: 
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Fi = mc^AcFc + (1) 

F2 = rnc, 2 Ac:>Fc + 111^^2 A^F^ ( 2 ) 

The only known quantities in these equations are the observed spectra, Fi and F 2 . Although 
we don’t know the intrinsic magnihcation factors of each component, it is possible to measure 
their ratios. In the absence of differential absorption or reddening these are just the wavelength- 
dependent flux ratios for each component. We set tjq = mcq/mc ^2 and rh = rriL^i/F the 
magnification ratios in continuum and emission line. We can then solve for C and C for one lensed 
image: 


Cl = mcqAcFc = ———(Fi - r]i^F2) 
Vc-Vl 


( 3 ) 


C-i — lA^Fi — ———(Fi — rjQF2) 

Vl-Vc 


( 4 ) 


We can think of these equations as representing a scaled subtraction of the spectrum of one 
lensed image (F 2 ) from another (Fi), with the scaling factor r} giving complete subtraction of one 
spectral component. In Eq. [3l for example, F 2 is scaled so that its emission line flux matches that 
of Fi. As long as the line-to-continuum ratio is different between Fi and F 2 then Fi — rji F 2 
eliminates the BEL spectrum leaving a representation of the absorbed continuum spectrum with 
a slope somewhere between that of lensed images 1 and 2. The scaling factor Tlc/iVc ~ Vl) 
then recovers the normalization and slope of the continuum spectrum seen in image 1. A similar 
interpretation of Eq. [J] eliminates the continuum component to give the lensed, absorbed BEL 
spectrum for image 1. 


Eor our decomposition of H1413+117 we define Fi to be the spectrum of image D and F 2 
to be that of image A. We take to be the wavelength-dependent continuum ratio from the 
fits described in Section [3l As seen in SectiSj the BEL magnification ratios are not wavelength 
dependent over each line and vary relatively little between Siiv, Civ, and Cm]. To reconstruct 
spectral components over the full wavelength range of the spectra we take t/l to be a constant: 
r]L = 0.40±0.02, the mean continuum-subtracted broad line flux ratio from Tabled] Eor individual 
BEL decompositions we use their individual magnification ratios (see below). 


Eigure [5] shows the resulting decomposition over the full spectrum. The upper panel shows 
Cd- the lensed, absorbed quasar continuum (Eq. [3]), with the continuum magnification factor of 
image D. We overplot the continuum fit for the image D spectrum. The lower panel shows Co- the 
lensed, absorbed BEL spectrum (Eq. d]), again as magnified in image D. We overplot the observed 
sp ectrum of image D minus the continuum fit and also an iron emission spectrum from the templates 
of IVestergaard fc WilkesI (1200 ll ). scaled to roughly match the iron features in the spectrum. 







- 17 - 


Figure [6] shows a nearly identical decomposition as Fig. [5] over the range of the Civ line, with 
the only difference being that we take rji = 0.410 ±0.003, the continuum-subtracted flux ratio over 
the red wing of the CiV line. Fig. [6] also shows the range of uncertainty for these decompositions 
(red lines): the maximum and minimum decomposed BEL and continuum components obtained by 
allowing the continuum flux to vary by ±5% (see Sect. [3]) and r]]^ to vary over its uncertainty range. 
It can be seen that the general shape of the decomposed components does not change signihcantly 
over the range of possible decompositions. 



Rest Wavelength (A) 


Fig. 6.— Upper. Absorbed continuum spectrum of lensed image D in H1413±117 around CiV {Cd 
from Eqn. [Sj black). Green dashed line shows the continuum fit to image D. Lower, absorbed Civ 
BEL spectrum of lensed image T) {Cd from Eqn. [H black). The green line shows the observed 
sp ectrum of image D with conti nuum fit subtracted. The blue line shows the iron emission template 
of IVestergaard fc WilkesI (j200l|) with arbitrary scaling. Red lines in both panels show the upper 
and lower limits of range of possible decompositions based on the uncertainti es in the continuum fit 
and in easured flux ratios. The blue line shows the iron emission template of IVestergaard fc Wilkes 
(j200l|) . Letters are designations for absorption trough features discussed in Section 15.31 


Lens image D is experiencing by far the strongest microlensing between the continuum source 
and BELR. However we do see some differences in broad line to continuum flux ratios in A, B and 
C. Figure [7] shows these same decompositions over the Siiv and Civ lines for image pairs A & D and 
for B & C. In these cases we take the A/D and B/C emission line flux ratios to be the average for 
the Siiv and CiV lines for each image pair. We see very consistent results, albeit with significantly 
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more noise in the B & C decomposition due to the smaller differences between these spectra. 
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Fig. 7.— Decomposed Siiv to Civ BEL {upper) and continuum (lower) spectral components of 
H1413+117 derived from lensed images A and D (black) and from lensed images B and C (red). 

If we divide Cd, the lensed, absorbed continuum spectrum, by our estimate of the lensed, 
unabsorbed continuum mc,D Fc we recover an estimate of the intrinsic absorption profile affecting 
the continuum. Figure [8] shows this for the Civ, Siiv, Nv, and Aim absorption troughs. 
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4.2. Comparison of Decomposition Approaches 


The derived spectral c omponents are s i milar to those found by HIO, and the derived absorption 
profile is similar to that of lAngonin et al.l (jlQQOl i. We discuss briefly similarities and differences in 
these works. 


HIO decompose in terms of macrolensed-only (their Fm) and microlensed + macrolensed (Fm^i) 
components, and use the combined spectra of lensed images A and B (Fab) as well as that of image 
D (Fjo). The spectral subtraction scaling factors used by HIO are M, the magnification ratio of 
the macrolensed-only component, and A, the magnification ratio of the microlensed -|- macrolensed 
component. A is taken as the continuum ratio in regions unaffected by broad line flux, and so 
is similar to our 'qQ. They measure a single value of A for each emission line from the average 
continuum bounding that line. This accounts for changes in the continuum magnification ratio 
from one line to the next, but not for variation across the wavelength range of the lines, however 
the latter is of order 1-2% and so is not a large effect. M is determined under the assumption that 
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at least part of each broad emission line is not microlensed, so that within the wavelength range 
of each line there is a region of the reconstructed micro+macrolensed component {Fm^) that is 
equivalent to the continuum spectrum. This can be interpreted as subtraction of Fab from Fb, 
with Fab scaled to eliminate at least part of the emission line (the part of the line closest to the 
continuum magnification). In practice, because the emission line magnifications are constant across 
their wavelength range (see Sect. ED, this should mean almost complete subtraction of emission 
lines and produces a result that is very similar to our Equation 2. 

(ADO) describe their decomposition in less detail than HIO, however 
their approach is mathematically similar to both ours and that of HIO. They express the spec¬ 
trum of lensed image D to be that of the (purportedly) unmicrolensed A-|-B-|-C composite plus 
a microlensed, absorbed continuum component. After deriving a continuum spectrum from the 
A-l-B-l-C composite they fit parameters to yield the absorption profile. This appears to be equiva¬ 
lent to performing a scaled subtraction of A-|-B-|-C from D to eliminate the BEL component and at 
the same time scaling out the wavelength-dependent continuum, reddening, lensing magnifications. 

In summary, although A90 and HIO formulate decompositions differently and make some differ¬ 
ent assumptions, all three approaches are mathematically similar and yield similar decompositions. 
It is worth noting that in our formulation the interpretation of the decomposition does not depend 
on their being any macrolensed-only component, only that the continuum-BELR ratio differs be¬ 
tween the two lensed images. While both HIO and ADO assume that the BELR is not microlensed, 
any such wavelength-independent microlensing or variability-related differences in the BELR do 
not invalidate the decompositions. 

We compare the details of these decomposition results in Section [5l 


Angonin et al.l (1990) 


5. Interpretation of Decomposition Components 
5.1. Interpretation of Broad Absorption Characteristics 


The continuum spectral component (Eig. [5]&:[6l upper panels) and the absorption profile that 
we derive from it (Eig. [HD provide a picture of the BAL absorption along the line-of-sight column 
to the accretion disk, and the shape of these profiles allows us to draw specific conclusions about 
the structure and kinematics of the outflow. We focus on interpretations in the context of rotating 
disk-wind models, in which the material responsible for both BAL and BEL f eatures arises from 
the accretion disk and is driven out by radiat ion and possibly gas pressure ([Murray Chiang 
19951 : IProga. Stone &: Kallmanll200Cll : lElvisll2000l i . In the case of Civ the continuum is experiencing 
essentially no absorption until ~1,500 km/s blueward of the emission line central wavelength, at 
which point transmission drops sharply over ~600 km/s to opacity ~ 1. This near-black region 
extends to ~7,500 km/s before transmission ramps roughly linearly to the end of the absorption 
feature at ~ 15,000 km/s blueward of line center. Two narrow emission peaks are seen in this latter 
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region. In the case of SiiV, absorption also has a sharp onset blueward of the Siiv AA1394 line 
center, in this case at ~3000 km/s. The deep region of its absorption trough ends much more 
sharply than in the case of Civ at ~7,500 km/s. A single clear emission peak is observed in the 
deep part of the trough. A weaker absorption trough appears redward of the deeper trough onset 
and may result from Siiv AA1403. In any case, all absorption is blueward of the Siiv doublet. As 
with Civ, the Nv trough begins at a blueshifted ~1,500 km/s, while the Aim trough appears at 
~4,500-5,000 km/s, though may show some absorption from ~2,000 km/s. 


Of particular interest is the sharp, blueshifted onset of this absorption, which is especially 
striking in Civ. This sharpness is not ambiguous, and is apparent over the full range of uncer¬ 
tainty in the decomposition (Fig. [6l lower). While BALs with highly blueshifted absorption onsets 
are sometimes observed (so-called detached troughs, blueshifted at several thousand km/s; e.g. 
Turnshek 1987), our decomposition of H1413-I-117 provides the first unambiguous description of a 
sudden, blueshifted absorption onset that is still within the velocity range of the BEL. The steep¬ 
ness and blueshift of this onset tells us two things: 1) the outflow has already accelerated to some 
degree before either it enters our sightline or it enters the region of efficient absorption within our 
sightline; and 2) the kinematics of absorbing gas in our sightline is strongly dominated by outflow, 
not orbital motion. 


Murray et al.l (j 19951 ) intepret highly detached absorption troughs as resulting from sightlines 


that graze the upper layer of an equitorial outflow, hence missing the yertica. 1 base of their wind 
model where the flow has a small radial component. iHigginbottom et al.l (j2ni3l ) (H13) apply radia¬ 
tive transfer code to their disk-wind model to predict Civ BAL profiles, also finding that these same 
grazing sightlines produce a very sharp, blueshifted onset. In these equitorial disk-winds we expect 
any deeper line of sight—those passing near the launch region of the wind—to show absorption at 
low velocities. Not only do they encounter launching gas prior to significant outward acceleration, 
but they also encounter gas with very high rotational velocities in the vicinity of the accretion disk. 
This latter effect is important, as shown by H13; at launch, the BAL flow shares the Keplerian 
velocity of the accretion disk and so is dominated by rotational motion. If the inner disk is occluded 
by this gas along our line of sight, and if this gas is significantly absorbing in the vicinity of the 
accretion disk, then some fraction of the continuum will be absorbed by the component of the 
orbiting flow that is receding along our sightline, producing partial absorption of the continuum 
at lower and/or redshifted projected velocities. On the other hand, even large rotational motion 
far from the accretion disk will only have a transverse velocity component within our line-of-sight 
column to th e disk. This argument against a rotational contribution to the line-of-sight velocity is 
supported bv iBorguet &: HutsemekersI (j2010|), who show that, for a spherically expanding outflow, a 
sharp absorption onset is only expected in cases with small rotational components. Any significant 
rotation results in a gradual absorption onset. 


This same grazing sightline in the H13 model also appears to reproduce the ramped blue edge 
of the Civ absorption trough that we observe in H1413-I-117. The models of iGiustini fc Proga 
(j2012l l show that a ramped increase in transmission at the blue edge indicates a large variation in 
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both velocity and density in the absorbing region along the line of sight, which appears to occur 
when this line of sight is close to tracing a streamline in the flow. This is consistent with a line of 
sight along the upper surface of the wind. 


The interpretation of iFiliz Ak et al.l (j2f)14l ) contradicts this. They classify SDSS BAL quasars 
according to the appearance of SiiV and Aim troughs alongside the ubiquitous Civ trough, and 
interpret those containing all three (H1413+117 does, making it a LoBAL) as resulting from lines 
of sight on the underside of the outflow, intersecting the base of the wind. 


The argument of the grazing sightline does not explain the different velocity shifts we see in the 
absorption onset for different line species, which should be the same if caused by our line of sight 
exiting the outflow. While ion species with the highest ionization potentials. Civ (64.5 eV) and Nv 
(97.9 eV), have similar ~1,500 km/s absorption onsets, the onset for Siiv (45.1 eV) is ~3,000 km/s 
and the deeper part of the Aim (29.4 eV) trough starts at ~4,500-5,000 km/s (though may show 
weaker absorption from ~2,000 km/s). This increase in absorption onset with diminishing ionization 
potential may suggest that for Siiv and Aim we are in fact seeing very sudden jumps in transition 
rate as ionizing radiation and/or gas density reach optimal levels as a function of radius and 
hence outflow veloci ty. For example, thi s may be the radius where it is first possible for “Locally 
Optimized Clouds” (iBaldwin et al.lll995l ) to form due to the diminishing ionization parameter. 


A combination of the above may be needed to explain the onset of all absorption troughs; we 
propose that for CiV and the higher ionization Nv we are seeing our line of sight exit the outflow 
near its launch point, while for lower ionization lines we are seeing the onset of favorable absorption 
conditions further out in the flow. In the ionization/radiative transfer calculations of H13, high Civ 
fractions are expected at the highest ionization parameters in their luminous benchmark quasar 
(Mbh = 10® Mq, Eddington ratio = 0.2), and so it is expected that Civ absorption should be seen 
at its inner-most radii. 


5.2. Interpretation of Emission Line Spectrum 


Eigures [5] & [6] (lower) shows a plausible BEL spectrum, however it should be remembered 
that this component includes the BAL absorption of this broad line emission. It is not possible to 
extricate the true absorption profile from the line profiles In the case of the Civ line (Fig. [6l 
lower) we see a region of near-zero emission in the rest frame 1515-1525A range (~4500~6400 km/s) 
where we might otherwise expect to see broad-line emission—either Civ or Fell, given the presence 
of one or both surrounding this reg ion, and the significant iron emission in this wavelength range 
of the IVestergaard fc WilkesI (120011 ) template. We take this as evidence of absorption of broad-line 
emission by the BAL wind. If this black region is due to absorption then we are seeing very high 
column densities and very high covering factors for absorption of the Civ and/or Fell emission at 
these velocities. This is very interesting due to the presumably much larger extent of the source of 
this broad-line emission compared to the continuum—the line-emitting gas is completely enveloped 
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by absorbing gas across the entire BELR at these projected velocities. 


Further evidence of absorption of BEL light by the BAL wind comes from the significant 
absorption of Lya by the Nv trough in this source (see also Turnshek 1987). While the symmetry of 
the Civ emission in this component resembles that of an unabsorbed BEL, it should be remembered 
that Civ very frequently shows a strong blue excess which may be absorbed in this case. The Siiv 
line shows a sharp drop in flux in its blue wing that also suggests absorption of BEL light. 


The iron emission template of IVestergaard fc WilkesI (120011 ) shown in these hgures provides a 
good match to several features in the emission line spectrum, and an imperfect match to others. 
The complexity and potential variation in iron emission makes a perfect match unlikely, however 
the scale and shape of these features match the template adequately. As well as natural variation, 
it is also possible that we are seeing differential microlensing acrosss the iron emission spectrum. 
The strong presence of iron features in component C implies that some of this emission arises on 
size scales similar to other BELs. 


The fact that the iron template provides a reasonable match to the iron emission across the full 
wavelength range of the BELR component (Fig. [5]) gives us conhdence in our choice of continuum, 
and that we avoided most of the iron emission; if not then we would expect a significant and 
wavelength-dependent mismatch between the decomposed BEL spectrum and the template. 


5.3. Features in the Civ and Siiv Troughs 

Low-velocity-width emission peaks are seen within the region of the Civ BAL absorption 
trough, marked in Figure (6) in the continuum component centered at blueshifts of ~8,000 km/s 
(a)) and ~11,000 km/s {h) and in the BEL component at ~4,000 km/s (c), ~7,500 km/s (d), and 
~13,000 km/s (e). An emission peak is also observed blueward of Siiv AA1393.76 at ~4,500 km/s 
(which we designate feature /, Fig. [8]). 

HIO observe a similar emission peak in their decomposition; it appears in all epochs in ~6,500— 
9,500 km/s blueward of Civ in their Ym (BEL) component (from the decomposition of lensed image 
D with respect to A &: B). This appears to correspond to our feature d. Also notable is a feature 
in their Fm^i (continuum) component from the decomposition of image C which corresponds to 
our feature a. A similar peak was also observed in the Civ absorption profile derived by A90 at 
~5,500 km/s. In all of these previous decompositions, signal quality makes it is difficult to rule out 
the presence of any of the other features that we note in the current work. 

The features observed in the Civ and Siiv troughs of the continuum component (a, b ^ f) 
indicate regions of reduced absorption in the BAL outflow. This may be due to radial fluctuations in 
overall outflow density and/or relative ion species density due to variations in ionization parameter. 
This interpretation reflects that of HIO, who propose a strongly absorbing cloud embedded within 
the outflow but with lower velocity width that results in a double-peaked Civ line. In the context 





- 23 - 


of a sightline grazing the upper lay er of a fast, equito r ial ou tflow, iGiustini &: Progal (j2012l i show 
that, in the hydrodynamic model of lProea &: KallmanI (j2004l i. “blobs” of failed outflow can cause 
dumpiness in this layer as they fall from the over-ionized high-inclination regions. 


Interpretation of the features in the BEL component is more complicated. Features d and e 
may result from clumpy absorption of a strong blue excess in Civ and/or Fell emission (which may 
or may not be absorbed). As discussed in Section [Ql the near-black region between features c and 
d does seem to indicate absorption of BEL light by the outflow. Whether this is also true of the 
dark region between features d and e is unclear, as feature e may be isolated Fell emission. Feature 
c is probably iron emission. 


It is clear that features a and d correspond to the same peak in the original spectrum, and 
so to the same region of diminished absorption. Feature b in the continuum component has no 
obvious counterpart in the BEL component. This may be due to intrinsically low BEL flux in this 
region or a difference in the BEL vs. continuum absorption profiles. 


5.4. Line-of-Sight Differences 


The differences observed in the BAL profiles between lens images cannot be due to differences 
in the line of sight. The BAL absorption profile affecting the continuum is determined by the 
region of the BAL flow sampled by the column along the line of sight to the accretion disk, and 
so different lensed images will only experience different absorption if their line-of-sight columns 
diverge signihcantly within this outflow. The image separation between images A and D is T.C8, 
which corresponds to a divergence angle of the A and D lines of sight from the source of C/.Tb and 
2'.'0 for zl = 0.94 and 1.4 respectively. Assuming zl = 1-4 and a conservatively small accretion 
disk with half-light radius 10^^ m (SS disk with L/Ledd = 0.1 and Mbh = 10® Mq), lines of sight 
have diverged by one half-light radius at ~100 light years. This is probably much larger than the 
strongly-absorbing region of the BAL flow. Choosin g a smaller zj gives an ev en smaller divergence. 
On the <1000 light day size scale of a large BELR (jKasDi et al.ll2000l . 120071 1. line-of-sight columns 
to the accretion disk still have > 98% overlap. 

Line-of-sight columns to the BELR, assuming 100 light day cross-section, diverge completely 
at ~28,000 LY, and so may be subject to different line-of-sight absorption features both within the 
host galaxy and in the ICM. This is also true of the continuum, however the continuum will be 
subject to differences in absorption features from a greater depth within the host galaxy. 

These line-of-sight di fferences do account for the many differences in narrow absorption features 
observed in our spectra. Monier. Turnshek fc Lup ie (jl999l ) show that these features arise both in 
the ICM and potentially the leasing galaxy, but probably not the host galaxy which has a higher 
redshift than observed features. We don’t investigate the absorption features in detail in this paper, 
however we note that the z = 1.43 Fel l lines {XreRt 2344. 2A, 2374 . 5A, 23 82. 8A) follow the same 
behavior in our spectra as observed by iMonier, Turnshek fc Lupid (jl999l L it is present in lensed 
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images A and B and absent in C and D. Several other features exhibit similar behavior, and we 
will investigate these in a separate study. 


6. Conclusions 

We have presented new spectroscopy of the gravitationally lensed broad absorption line quasar 
H1413+117 using Gemini North’s CMOS IFU. Extracted spectra of the four lensed images reveal 
that differential microlensing strongly affects the relative magnification of the broad emission lines 
compared to the continuum in lensed image D, with the continuum in D appearing anomalously 
bright compared to all other lensed images. Lensed images A, B, and C also show small differences 
in their continuum-to-broad line flux ratios relative to each other. 

Subtracting a fit to the continuum emission outside the broad emission and absorption regions 
allows us to determine the wavelength-dependent D/A flux ratios of the broad emission lines. These 
flux ratios are flat, indicating no differential microlensing across the velocity structure of the BELR. 

Monte Carlo simulations were performed to analyze the dependence of the flatness of broad 
line flux ratios on the the size range spanned by the BELR velocity structure. It was found 
that BELRs with an upper size limit in the 10^^—10^^ m range, and spanning 0.5 dex across their 
velocity structure are very unlikely to produce the observed flat flux ratios. Smaller BELRs are also 
unlikely due to the strong differential magnification observed between broad lines and continuum. 
Therefore, it is inferred that the BELR in H1413-I-117 is significantly larger than 10^^ m, and 
produces no significant emission from regions smaller than ~ 5 x 10^^ m, or 20 light days. 

We take advantage of the differential magnification between lensed images A and D to alge¬ 
braically decompose the quasar spectrum into two components: the absorbed broad line spectrum 
and the absorbed continuum spectrum, and we use the latter to derive intrinsic absorption profiles. 
Our decomposition leads us to these conclusions; 

• A sharp, blue-shifted onset of absorption is observed in all prominent BAL features in the con¬ 
tinuum component. In the case of Civ, Nv, and Lya this onset is at ~1,500 km/s, which may 
represent the inner edge of a radial outflow. In the context of typical disk-wind models this would 
place the onset of BAL absorption near its launch point, and so in the vicinity of the accretion disk. 
In the case of Siiv and Aim, which have a larger absorption onsets, we may be seeing increases in 
absorption efficiency of lower ionization species with distance from the ionizing source. This sharp 
onset also indicates that the line-of-sight velocity is dominated by outflow, not rotational motion. 

• Overall the derived absorption profile for Civ resembles that of an accelerating outflow that is 
either spherical or viewed relatively close to a streamline, albeit with a ~l,500km/s initial velocity. 

• The Civ and/or Fell BELR appears to be strongly absorbed in the range ~4,500—6,400 km/s. It 
is unclear whether emission features blueward of this trough are Civ or Fell, however the complete 
absorption of either of these features suggests a very high covering factor for the absorbing material 
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at large radii. Combined with the apparent strong absorption near the outflow launch point, this 
suggests that gas responsible for BAL features maintains a high-covering-factor over at least two 
orders of magnitude in size. 

• We observe regions of reduced absorption within the BAL outflow, indicating dumpiness in gas 
density or relative ion species density. 
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Fig. 8.— Intrinsic absorption spectrum of BAL wind for (from top to bottom) Civ, Siiv, Nv, and 
Aim based on continuum spectrum derived from Eqn. [3l Velocity zero-point is set at line center 
or at lower wavelength line center for doublets. Vertical dotted lines show line centers. Letters are 
designations for absorption trough features discussed in Section 15.31 





































































































































































